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Abstract. We present adaptive optics assisted OASIS integral field spectrography of the SO galaxy NGC 4621. Two- 
dimensional stellar kinematical maps (mean velocity and dispersion) reveal the presence of a ~ 60 pc diameter counter-rotating 
core (CRC), the smallest observed to date. The OASIS data also suggests that the kinematic center of the CRC is slightly offset 
from the center of the outer isophotes. This seems to be confirmed by archival HST/STIS data. We also present the HST/WFPC2 
V - I colour map, which exhibits a central elongated red structure, also slightly off-centered in the same direction as the kine- 
matic centre. We then construct an axisymmetric mod el of N GC 4621: the two-integral distribution function is derived using 
the Multi-Gaussian Expansion and the Hunter & Qian (1993) formalisms. Although the stellar velocities are reasonably fitted. 



including the region of the counter-rotating core, significant discrepancies between the model and the observations demonstrate 
the need for a more general model (e.g. a three-integral model). 
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1. Introduction 

The fact that early-type galaxies are not necessarily oblate 
major-axis rotators has been well known since the end of 
the 1970s. The number of objects showing signifi cant m inor- 
axis rotation attained ten in 1988 (Wagner et al. 



1988), and 



some of these objects later showed evidence for kinemat- 
ically d ecoup led core structures (NGC 4365, NGC 4406 : 
Bender |l988l NGC 4589: Moellenhoff^ & Bender |l989; 



NGC 5982: Wagner 1990). The most cited formation scenario 
for such structures is the hierarchical model, which involves 



single or multiple merging steps (see e.g. Cole et al. 1994 and 
references therein). Observations of kinematically decoupled 
cores can therefore be used to constrain the merger tree of 
galaxies. However the total fraction of early-type galaxies con- 



taining such a core is poorly known (de Zeeuw & Franx 1991 ) 



Such cores are generally revealed using long-slit measurements 
with 1" FWHM or worse spatial resolution: this introduces 
an a priori assumption on the geometry of the central struc- 
tures, and prevents the detection of cores with small apparent 
sizes. Finally, the rough characteristics of these cores (fraction 
in mass, relative age and metallicity) remain very uncertain, as 
very few studies exist and are often limited to morphological 



data (see Carollo et al. 1997 for 15 KDCs) 



Send offprint requests to: F. Wemli, e-mail: 
fabien@obs .univ-lyonl . fr 

* Based on observations collected at Canada-France-Hawaii 
Telescope, which is operated by CNRS of France, NRC of Canada, 
and the University of Hawaii. 



Recent studies (e.g. Verolme et al. 2002 ) show that state-of- 
the-art models combined with two-dimensional integral-field 
spectroscopic data are required to constrain precisely the global 
physical parameters of early-type galaxies, such as the inclina- 
tion angle or the mass-to-light ratio, as well as the central char- 
acteristics for objects exhibiting features such as kinematically 
decoupled cores (e.g. the core mass of IC 1459, Cappellari et 
al. 2002). Moreover, integral-field spectroscopy also provides 
access to line-strength maps: the detailed kinematic informa- 
tion can thus be coupled to the line indices in order to more 
easily disentangle the core from the host galaxy (Davies et al. 
2001). 



The SAURON survey (Bacon et al. 2001a, de Zeeuw et 
al. 2002) was carried out in order to study the integral field 
kinematics of a sample of 72 early-type galaxies over a wide 
field of view (> 33"x41"). In parallel, to obtain complemen- 
tary high spatial resolution data, we designed a program to 
observe their central parts (~5"x5") with the integral field 



spectrograph OASIS (Bacon et al. 2000). The aim is to probe 



the central spatial morphology and dynamics of a subset of 
galaxies in the SAURON sample, in order to link them with the 
global properties of the host galaxies provided by SAURON and 
oth er wid e field studies. The additional use of PUEO (Rigaut et 
al. 1998), the Canada France Hawaii Telescope AO bonnette, 
with OASIS allows us to reach resolutions of about 0'.'25 
(~ 12 pc at 10 Mpc) in the red part of the visible when bright 
guiding sources are available. 



NGC 4621, an SO galaxy (Lauer et al. 1995), with My = 



9.6, located in the Virgo cluster (D=18.3 Mpc, Tonry et 
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al. 2001 1, is the first target of the subsample that was observed 
with OASIS/PUEO. Its steep power-law central luminosity pro- 
file (y - 2.03, Gebhardt et al. 1996) is ideally adapted to serve 
as a guiding source for the AO. NGC 4621 is close to edge- 
on, and its innermost 6" isodensities reveal a nuclear disk 
(Sil'Chenko et al. 



1997) 



In this paper, we present the first results of our OASIS ob- 
servations of NGC 4621, with the discovery of the smallest 
counter-rotating stellar core observed to date. The scale of the 
core is 60 pc, while the previously detected cores have an av- 
erage radius of around 1 kpc (see the sample of Carollo et 
al. 1997 ). Detailed information on the different data sets used in 



this work, including HST/WFPC2 and STIS observations, are 
provided in Sect. ^ The corresponding measured photometry 
and stellar kinematics are presented in Sect. ||. We computed 
a two-integral distribution function (DF) model of the galaxy, 
which is presented in Sect. ^ The discussion is carried out in 
Sect. 15. 



2. Observational data sets 

2.1. OASIS data 
2.1 .1 . Observations 

Ca-triplet OASIS exposures of NGC 4621 were obtained in 
January 2000 using the medium spectral configuration (MR3: 
8346-9152 A) at the f/20 (PUEO) focus of the CFHT, and a 
spatial sampling of 0'.'16. Eight out of ten 30 mn exposures 
were fully reduced using the dedicated XOasis package]^ and 
merged into a single frame (two frames were discarded due to 
bad seeing conditions and associated loss of guiding). The re- 
sulting 3D dataset0 probes about 4" x 4" (300x300 pc^), with 
a resolution of cr = 70 km.s '. We modeled the PSF of the 
merged datacube using the sum of two concentric Gaussians: 
their parameters were retrieved (via a least-square procedure) 
by comparis on wit h the WFPC2/HST F814W frame as in 
Bacon et al. ( |2001b[ ). The resulting FWHM of the merged dat- 
acube is Off 5 1 . 

The total S/N (summed over the bandwidth) drops very 
quickly towards the outer parts, as it is 65 at the center, 20 at 1" 
and 10 at 2"on the major-axis. We thus convolved the datacube 
with a Gaussian of variable width in order to obtain a higher 
S/N over the field (see Fig. |l]). 



2.1.2. Kinematic measurements 

We used a tuned version of the Fourier Correlation Quotient 



(FCQ; Bender 1990 1 method to derive the line of sight veloc- 



ity distributions (LOSVD) at each measured spatial position. 
We used a single template star (G8III: HD073665), but tested 
the stability of the kinematical mesurements with different tem- 
plates, and different continuum subtraction parameters, only 
resulting in minor differences. The LOSVDs were then fitted 
using a single Gaussian function, yielding the mean velocity 
and dispersion maps. We could not obtain reasonable maps of 



higher order moments, due to the rather low signal to noise ra- 
tio (even after the convolution performed on the datacube). The 
error bars on the velocity and dispersion measurements were 
computed via a Monte-Carlo algorithm using the measured S/N 
(see Copin|200C|, PhD Thesis). 



2.2. Ancillary data 
2.2.1. Pliotometry 

In order to obtain an accurate three-dimensional mass model 
of the galaxy we need both wide-field and high-resolution 
photometry to probe the visible mass up to large radii and 
to correctly sample the central structure (the cusp). We used 
a 20' V band image of the galaxy taken at the OHP 2m 
Telescope, kindly p rovide d by R. Michard (Idiart, Michard, & 
de Freitas Pacheco 2002). The central regions were examined 
with the help of WFPC2/HST data^ retrieved from the archive 



(F555W and F814W filters, Faber et al. 1997, #5512). For 



both bands, three unsaturated frames were merged and cosmic 
ray-corrected (330 s and 230 s exposure time for F555W and 
F814W filters respectively). We performed photometric cali- 
bration using the VEGAMAG standard, and used the PSFs com- 
puted using tinytiiti^ We adjusted the levels of the OHP im- 
age to the F555W frame taking into account the different spec- 
tral domains and the PSF. 

2.2.2. Long-slit kinematics 



We used long slit data from Bender, Saglia, & Gerhard (1994, 
hereafter BSG94; data kindly provided by R. Bender) in or- 
der to complement the OASIS data up to large radii: major and 
minor axis velocity and dispersion measurements are available 
inside 40" (3.5 kpc), with a seeing of ~1'.'8 FWHM (other kine- 
matical - and line-strength - measurements were published by 
Sil'chenko 1997 but with a lower spatial resolution). 

We finally reduced unpublished HST/STIS major-axis 
data (Green, ID #8018; G750M grating, Ca-Triplet) in or- 
der to examine the central kinematics of NGC 4621 at high 
spatial resolution. Two exposures with a total of 72 mn 
were available. The calibrated data were retrieved via the 
STScI data archival system (calstis pipeline). Appropriate 
flat field exposures were retrieved to correct for the fring- 
ing (mkfri ngef l at/defringe IRAF routines, Goudfrooij & 
Christensen 1998), critical at these wavelengths, as well as five 
KO-III kinematical template star exposures (HR7615, Green, 
#7566). Further rejection of cosmic rays was then performed 
on the NGC 462 1 individual exposures, before recentring and 
merging. The stellar kinematics were extracted using the same 
FCQ and fitting methods as for the OASIS data. As for the 
OASIS datacubes, the S/N of the STIS data was found sufficient 
only to derive the mean stellar velocity and dispersion within 
the central arcsecond. 



' |ittp://www-obs.univ-lyonl . fr/" oasis/home/index .html 
^ (a, 6, A), one spectrum per spatial element 



based on observations collected with the NASA/ESA HST, ob- 
tained at STScI, which is operated by AURA, Inc, under NASA con- 
tract NAS5-26555 



* fittp : //www . stsci . edu/ sof tware/tinytim/tinytim . html 
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3. Results 



3.2. HST photometry 



To avoid any confusion due to the complexity in both the pho- 
tometry and the kinematics of NGC 4621, we will use the fol- 
lowing convention: the photometric major axis of the galaxy (as 
measured by HST) defines the x-axis, the negative values be- 
ing in the SE quadrant, referred to as the left side of the galaxy. 
The center (0, 0) is defined as the center of the HST isophotes of 
NGC 4621 within 1", excluding the central U.'l, the isophotes 
of which exhibit a significant asymmetry (see next Section). 
The centering and rotation procedures of the HST frames have 
been performed using an algorithm which minimizes the stan- 
dard deviation of the recentered and rotated frame subtracted 
by its flipped counterpart. This leads to accuracies of 0. 1 de- 
grees and (X'OOS. All figures in this paper share a common ori- 
entation, as shown in Fig. |[ 



3.1. Stellar velocity and velocity dispersion 

The OASIS stellar mean velocity and dispersion maps are pre- 
sented in Fig. ^ The datacubes have been centered and rotated 
in order to match the HST data, using the center and angle pro- 
vided by the PSF fitting procedure mentioned in Section^. The 
velocity field reveals a clear counter-rotating core (CRC). The 
position of the zero-velocity curve is used to measure the cen- 
ter and the size^of the CRC: it is I'.'l (150 pc) in diameter, and 
off-centered by ~ 0'.'2 towards the SE. The total velocity am- 
plitude observed within the central 1" along the major axis is 
100 km.s while the peak-to-peak velocity amplitude of the 
CRC's is only 35 km.s The stellar velocity dispersion map 
peaks at ~ 330 + 25 km.s"'. The dispersion map however ex- 
hibits high frequency substructures in the central 0'.'5 (Fig.|l]), 
the minor-axis dispersion profile even having a local minimum 
at the centre. Considering the amplitude of these structures and 
their spatial scales (comparable to the local FWHM), we should 
wait for data with better signal to noise ratio to discuss these 
features. 

The BSG94 data show no significant minor-axis rotation, 
and a maximum velocity of 140 km.s"' at 30"(see BSG94 or 
Fig. I?]). Besides a relatively weak velocity amplitude in the in- 
nermost 4" (< 60 km.s"'), these data do not show any hint 
for the existence of the counter-rotating core. We will use the 
BSG94 kinematics to constrain the dynamical model outside 
the OASIS field of view. 

The STIS velocities and dispersions are shown in Fig. ^ In 
these data the size of the CRC is Off! (~ 60 pc) and its veloc- 
ity center is located at jc = -0'.'05 (~ 4.5 pc). We have care- 
fully checked the centering of the STIS data with respect to the 
reference WFPC2 images, and found that this offset is robust. 
At STIS resolution the major-axis velocities of the CRC reach 
~ 75 km.s"'with respect to the systemic. The maximum veloc- 
ity dispersion is 320 + 27 km.s"', at x = 0'.'05. These values are 
consistent with the OASIS data considering the higher spatial 
resolution of the STIS data. 



Both F555W and F814W images reveal a peculiar structure in 
the core. Indeed, the luminosity peak is offset from the center 
of the outer isophotes by O'.'Ol towards the east. The photo- 
metric peak is located in the upper left quadrant (Fig. ||). This 
feature does not seem to be an artefact, as it is clearly visible 
on each individual HST frame, and in both bands, and cannot 
be attributed to centering uncertainties (0'.'005). 




< 
o 



I 



-0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 
arcsec 

Fig. 3. Inner isophotes of NGC 4621 (WFPC2 F555W, solid 
contours, step of 0.2 mag). The center (0,0) of the galaxy is 
defined as the center of the outer isophotes (outside O'.'l). The 
open triangle and filled circle correspond to the center of the 
counter-rotating core (CRC) as measured in the OASIS and 
STIS data repectively (see Sect. 



3.1). The vertical bars cor- 



The limits of the CRC are defined by the zero- velocity curve 



respond to the OASIS and STIS spatial resolutions (cr). The 
dashed contours represent the V - I colour map from WFPC2: 
levels 1.36 and 1.37 mag. North is 73 degrees clockwise from 
the vertical axis. 



The two HST frames were then PSF-crossconvolved and 
divided to obtain a V - I colour image (Fig. It reveals a 
central gradient (V-I increases towards the center) with a y- 
axis elongated structure centered at (-0'.'02, O'.'Ol) (see dashed 
line on Fig. ||). The isochromes have slightly higher ellipticities 
than the isophotes (s = 0.5 at 2'.'5 for V-I, vs. 0.4 for the V 
frame). 

4. Models 

In this Section, we present the methods used to model the pho- 
tometry and kinematics of the central region of NGC 462 1 . The 
available data suggest that the very central region (100 pc) of 
NGC 4621 slightly departs from axisymmetry. The CRC seems 
off-centered in both the OASIS and STIS data. However, the 
off-centering is only ~ 4 pc. At the resolution of the BSG94 
data, this is obviously not resolved. Moreover, the position an- 
gle measured in the OHP photometry does not vary more than 2 
degrees within 15'. Axisymmetry is therefore a reasonable ap- 
proximation for a first modelling. Axisymmetric, two-integral 
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Fig. 1. OASIS stellar kinematic data of NGC 4621. Top left: mean velocity field. The zero-velocity curve has been emphasized 
with a white thick line. Top right: velocity dispersion. White contours are reconstructed isophotes. The data have been smoothed 
to improve the signal to noise: the black circles represent the beam sizes (FWHM). Below the cuts along the major and the minor 
axes are shown along with the best fitting two-integral model (dashed lines, see Sect. 0). 
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Fig. 2. STIS mean velocity (left) and dispersion (right) profiles along the major axis. The dashed line corresponds to the best 
fitting two-integral model (see Sect. The velocity dispersion is clearly underestimated by the present model, which hints for 
the need of an additional central dark mass, a more general treatment (three-integral model) or both. 





1.27 




Fig. 4. V - I map (WFPC2, F555W and F814W filters). Left: convolved with a gaussian of cr = (7.' 1. Right: convolved with the 
PSF of the OASIS data ((X'Sl FWHM). Contour separation is 0.005 mag. A y-axis elongated central structure can be observed at 
HST resolution (see also Fig. ||). 
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Fig. 5. Top panels: Multi Gaussian Expansion fit (thick contours) of NGC 4621 superimposed on the V band isophotes (thin 
contours). Top left : OHP V-band photometry. Top right : HST/WFPC2 F555W band image (isophote step of 0.4 mag/arcsec^). 
Notice the nuclear disc in the HST data. Bottom left panel: NGC 4621 light profiles along - x^/a^ + y^/b^, where a and b are 
the semi major and minor axes of fitt ed ellipse respectively. Crosses correspond to the WFPC2 image, circles to the deconvolved 
WFPCl data from Byun et al. ( 1996 ). The dashed and dotted curves correspond to the convolved and deconvolved MGE models 
respectively. Bottom right panel: ellipticity profile, with Michard's (solid bold line) and WFPC2 (solid hairline) data, along with 
the MGE model (dotted line). 



models then have the advantage to be semi-analytical. We first 
derived simple Jeans models to roughly constrain the input pa- 
rameters (inclination, mass-to-light ratio). We then used the 



Hunter & Qian (1993) formalism to compute the two-integral 



distribution function / of NGC 462 1 , as a function of energy 
(E) and of the vertical component of angular momentum (L,). 
The present modelling is only intended to provide a first view 
at the dynamics of NGC 4621, so we decided not to include a 
central dark component: this issue will be examined in detail in 
a forthcoming paper 



4.1. MGE 

The first step in modelling NGC 4621 was to build a lumi- 
nosity model which properly reproduced the observed pho- 
tometry. We used the Multi Gaussia n Exp ansion (MGE) for- 
mal ism p roposed by Monnet et al. (1992) and Emsellem et 
al. ( 1994 ), which expresses the surface brightness as a sum of 
two-dimensional Gaussians. Assuming the spatial luminosity is 
also a sum of (three-dimensional) Gaussians, given the choice 
of viewing angles, and using an MGE model for the PSF, we 
could then deconvolve and deproject the MGE model uniquely 
and analytically. 

The procedure was performed stepwise. We began fitting 
the wide field V-band image. We then subtracted the outer 



F. Wernli et al.: A 60 pc counter-rotating core in NGC 4621 



7 





liLQ.pc---) 




g 


1 


4.502 10^ 


0.040 


0.860 


2 


7.713 lO-* 


0.112 


0.610 


3 


4.792 lO'* 


0.201 


0.941 


4 


1.161 lO'* 


0.438 


0.344 


5 


1.475 lO'* 


0.516 


0.919 


6 


4.638 10^ 


1.036 


0.325 


7 


6.317 10^ 


1.280 


0.872 


8 


3.357 10^ 


2.486 


0.275 


9 


2.700 lO"* 


3.211 


0.658 


10 


1.326 lO"* 
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12.468 
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14 
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0.633 


15 


1.517 10' 
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a=inf 



Table 1. MGE photometric model of NGC 4621 in the V band. 
The flattened inner disc components are emphasized in bold. 



gaussian components from the high resolution WFPC2 image, 
and fitted the residual image (central 15")- In this step, we had 
to exclude the innermost 0'.'4, to avoid convergence problems 
probably due to the slightly asymmetric central feature (Fig. |]). 
We finally fitted the very central arcsecond. Gathering these 
three parts, we thus obtained a 15 gaussian components model, 
with the same center and PA, the parameters of which are given 
in Table |lj The goodness of the fit is illustrated in Fig. || 

4.2. Two-integral models 



We then made use of the Hunter & Qian (1993) formalism to 
build the two-integral distribution function of the galaxy using 
the best fit value for the mass-to-light ratio of Ty - 6.6 Mq/Lq 
found from simple Jeans models. We used a default value of 
/ = 90°, which produced a marginally better fit. 

The DF is divided into two parts, which are respectively 
even and odd in L,. The even part (/even) is uniquely determined 
by the input MGE mass model. This involves the calculation 
of a path integral in the complex plane, as described by Hunter 
& Qian (1993). The odd part of the distribution function (/odd) 



is then parametrized, and adjusted to fit the kinematical data. 
We chose the parametrization proposed by van der Marel et 



al. (1994), and modified it to account for the CRC. The origi- 
nal parametrization corresponds to Eq. [l|, and we additionally 
allowed j] to be function of Ep = E/E^^,^ (Fig. ||). The analyti- 
cal form of this t] function is the same as the one in Eq. ||, with 
an additional variable change. We can adjust the energy Ea,t 
above which the stars begin to counter-rotate, as well as the 
smoothness of the transition (a = oo: abrupt transition, a = 1: 
smooth transition, see Fig. o). 



foUE, = (2/7 - 1) /l„(LJLj,„„) ■ /even(£', L,) (1) 

( tanh(Q'jt:/2)/ tanh(Q'/2) (a > 0) 
ha(x) = i X (a = 0) (2) 

[ (2/a)arctanh[xtanh(a/2)] (a < 0) 



r\(Ep) 




Ep 

Fig. 6. Parametrization of rj. Ep is the normalized energy Ep = 
E/E„y^^. Eait is the energy at which the counter-rotation starts, 
e.g. when j] < 0.5, /odd < (see Eq. [l|). The bold line corre- 
sponds to the best-fitting model. 



The best fit model reproduces the BSG94 velocity profiles 
reasonably well (Fig. ^ with values of a ranging from 8 out- 
side 10" to -2 in the central part. The higher resolution OASIS 
(Figs. ^ and |l]), and STIS (Fig. ^ velocity measurements, re- 
vealing the counter-rotating core, are also well fitted by this 
two-integral model. The best fitting model uses a core with a 
diameter of 1" 1 (£„« - 0.62), and an abrupt transition (a = 100 
i.e. almost all stars having Ep > E^^, are counter-rotating, see 
bold line in Fig. A rough estimate of the CRC mass can 
be made by selecting stars counter-rotating in the central part. 
This is performed by integrating the DF weighted by a function 
which is for E < Ea-,t and > and 1 otherwise. The total 
mass of NGC 4621, which is given by the mass-to-light ratio 
and the deprojected MGE-model is 1.78 ■ 10" Mq. The total 
mass of the CRC is 2.13 • 10** Mq, yielding a mass fraction of 
0.12%. 

The dispersion profiles are well reproduced by the model 
outside the central few arcseconds. The central values of the 
dispersion predicted by the models are however systematically 
too low compared to the BSG94 observations. This is con- 
firmed by the OASIS and STIS dispersion values: the model 
predicts a central dispersion of ~ 220 km.s"', to be compared 
with the central observed STIS dispersion of 320 + 27 km.s"'. 
We were finally unable to fit the higher order moments, even 
at large radii. The values predicted by the model are thus 
systematically too high, by a factor of almost two. This dis- 
crepancy could not be solved even by changing the parameters 
of the odd part of the DF. 

These two discrepancies do indicate that we need a more 
general dynamical model for NGC 462 1 . First, we should re- 
move the constraint imposed by the two-integral model by al- 
lowing a third integral of motion. There may then still be the 
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Fig. 7. Best fit two-integral DF model. Left: Major axis. Right: 



Fig. 8. Best fit two-integral DF model for the comparison with 
OASIS. Only the velocity map is shown: central dispersion val- 



Minor axis. Inclination is 90 degrees, no black hole and Ty = ^^s are too low due to the absense of a central dai-k mass. Grey 
6.06. Points are the BSG94 data, and the curve is the model, levels range from -50 to -1-50 km.s ' 



The fit is reasonable but indicates the need for a black hole 
(low central cr) and maybe a third integral (overestimated hs). 



need for an additional central dark mass to explain the observed 
dispersion values. Such a model will be examined in a forth- 
coming paper (Wernli et al., in preparation). 

5. Discussion and conclusion 

In this paper, we report the discovery of a ~ 60 pc diame- 
ter counter-rotating core using new OASIS integral field spec- 
troscopy, a result confirmed with archival STIS long-slit data. 
A structure at this scale has never been observed before, except 
in M 31, whi ch exh ibits a 20 pc co-rotating decoupled core 
(Bacon et al. ZOO lb and references therein). Typical sizes of 
kinematically decoupled co res ar e at the order of 1 kpc (See e.g. 
the sample of Carollo et al. 1997 ). The computed mass fraction 
of 0.12% (2.1 • 10^ Mq) is on the low side, as compared for 
instance to the CRC of IC 1459, which is 0.5% for a core size 



center of the outer isophotes. This off'set seems consistent with 
the weak asymmetry detected in the WFPC2/HST V and / 
frames and the V - I colour map (Fig. ||). The central struc- 
ture elongated along the minor axis detected in the V-I colour 
map represents an increase of V - / = 0.03 which could either 
be due to an intrinsic stellar population gradient, or to dust. In 
the latter case, this would correspond to Av = 0.06 (assuming 
Ry - 3.1 and R j - 1.86 from models of the Galaxy, Rieke & 



of 3 kpc (Cappellari et al. 2002), but for which the estimation 
is of higher accuracy. As a comparison, the black hole mass 
from the cr - M. re lation for IC 1459 is 1.1 ■ 10'' Mq (Merritt 
& Ferrarese pOOTb| ) and 2.13 • lO'^ Mq for NGC 4621 (Merritt 
& Ferrarese 2001). The mass fraction estimate provided here 
in the case of NGC 4621 can only be indicative: the DF sepa- 
ration is inaccurate and difficult, and to get a precise result we 
would need an analysis in integral space. 

Both individual OASIS and STIS data sets suggest that the 
CRC of NGC 4621 is slightly off-centered with respect to the 



Lebofsky 1985 ). Note that we do not detect any high frequency 
structure in the individual WFPC2 V and / frames. 

The two-integral models reasonably reproduce the ob- 
served velocity profiles in the outer parts, as well as in the 
region of the CRC (OASIS and STIS). These axisymmetric 
models do obviously not take into account the observed off- 
centering. The predicted central value of the velocity dispersion 
is significantly too low by 50 km.s We were also unable to 
simultaneously fit higher order Gauss-Hermite moments: e.g. 
the predicted h^, is systematically too high (see Fig. This 
demonstrates that we need more general dynamical models in 
order to properly fit the observed kinematics. We are there- 
fore in the process of constructing three-integral Schwarzschild 
models, with the possibility to include a central dark mass. This 
may solve the discrepancies mentioned above. 

If one assumes a merger-scenario for the CRC's origin, the 
off-centering, if confirmed, could be either due to the fact that 
the merging process is still on-going, or may be the result of a 
stable mode (see e.g. M 31 in Bacon et al. 2001t). The short 
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dynamical timescale at the radius of the CRC (1 Myr) seems to 
favour the latter hypothesis. Further discussions regarding this 
issue must however wait for additional spectroscopy. A detailed 
study at high spatial resolution of the stellar populations in the 
central arcsecond would certainly help in this context. Existing 
high resolution STIS spectra at bluer wavelengths can be used 
for this purpose. 

Stellar kinematics at HST resolution is available for only 
a handful of early-type galaxies. M 31 is one of the rare ex- 
amples where it is possible to have sufficient spatial resolution 
to measure kinematical features at parsec scales, such as its 



20 pc nucleus (Bacon et al. 2001b and references therein). The 



OASIS data presented in this paper demonstrate that it is possi- 
ble to study substructures with a characteristic size of 60 pc in 
galaxies at the distance of the Virgo cluster The presence of the 
CRC in NGC 4621 raises new questions about the dynamical 
status of the centers of early-type galaxies. What is the frac- 
tion of early-type galaxies that have such substructures? Two- 
dimensional spectroscopy at high spatial resolution is clearly 
needed to simultaneously study the dynamical and chemical 
contents of these cores, and to eventually understand their ori- 
gin. 
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